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SIGNATURES OF EXTRAGALACTIC DUST IN PRE-SWIFT 

GRB AFTERGLOWS 

D. A. Kann,^ S. Klose,^ A. ZeU 

ABSTRACT 

We present the results of a systematic analysis of gamma-ray burst afterglow 
spectral energy distributions (SEDs) in the optical/near-infrared bands. Our 
input list includes the entire world sample of afterglows observed in the pre- 
Swift era by the end of 2004 that have sufhcient publicly available data. We 
apply various dust extinction models to fit the observed SEDs (Milky Way, Large 
Magellanic Cloud and Small Magellanic Cloud) and derive the corresponding 
intrinsic extinction in the GRB host galaxies and the intrinsic spectral slopes of 
the afterglows. We then use these results to explore the parameter space of the 
power-law index of the electron distribution function and to derive the absolute 
magnitudes of the unextinguished afterglows. 

Subject headings: dust, extinction — galaxies: high-redshift — gamma rays: 
bursts 


1. INTRODUCTION 

According to the most popular progenitor model of long-duration Gamma-Ray Bursts 
(GRBs), the collapsar model (Woosley 1993), a GRB is the result of ultra-relativistic jets 
ejected by an accreting black hole formed by the core-collapse of a massive star (most prob¬ 
ably a Wolf-Rayet star). This predicts a physical link between GRBs and supernovae (SNe) 
that has been spectroscopically confirmed in four cases so far: XRE 020903 (Soderberg et 
al. 2005), GRB 021211/SN20021t (Della Valle et al. 2003), GRB 030329/SN2003dh (e.g., 
Hjorth et al. 2003a; Stanek et al. 2003) and GRB 031203/SN20031w (e.g., Malesani et al. 
2004, Mazzali et al. 2005, in preparation). Eurther evidence comes from the observation 
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of late-time bumps in GRB afterglows that can be modelled very well by an underlying 
SN component (Bloom et al. 1999; Zeh, Klose, & Hartmann 2004), and which have led to 
the conclusion that in fact all long-duration GRBs are physically related to SN explosions 
(Zeh, Klose, & Hartmann 2004, Paper I). Furthermore, the collapsar model implies that the 
progenitors of long-duration GRBs are associated with regions of high-mass star formation 
(Paczyhski 1998), which might reveal themselves by a detectable extinction in the GRB host 
galaxies along the lines of sight towards the burster. This idea is further supported by the 
so-called dark bursts (Groot et al. 1998), for which no optical afterglow has been discovered 
despite rapid and deep searches in small error box regions. While most none-discoveries are 
the result of too shallow searches and too large error boxes (e.g., Jakobsson et al. 2004b; Rol 
et al. 2005), a small percentile remains that require intrinsic extinction to dim the afterglow, 
e.g., GRB 970828 (Groot et al. 1998; Djorgovski et al. 2001), 990506 (Taylor et al. 2000) and 
020819 (Jakobsson et al. 2005a), while others might have been intrinsically underluminous 
(for a discussion, see e.g., Fynbo et al. 2001a; Lazzati, Govino &: Ghisellini 2002; Klose et 
al. 2003). 

Of particular interest within the dust extinction model is the statistical distribution 
of the amount of visual extinction in the GRB host galaxies, as well as the nature of the 
corresponding dust. Given the fact that GRBs and their afterglows can be observed up to 
high redshift, they offer the possibility to get insight into the nature of cosmic dust when 
the universe was much younger and galaxies were much less evolved. Furthermore, since the 
optical properties of the dust grains trace the environmental conditions in the interstellar 
medium (Fitzpatrick & Massa 1986; Draine 2000; Bradley et al. 2005), they are to some 
degree an indicator for the physical conditions to make a GRB progenitor. 

The present paper is the third in a series of papers where we employ a large database 
of photometry gathered from the literature to reanalyze all optical afterglow light curves of 
GRBs in the pie-Swift era in a consistent manner to derive a homogenous sample, which is 
then used to study afterglow properties in a statistical way. In Paper I and an update (Zeh, 
Klose, & Hartmann 2005), the properties of the supernovae underlying nearby afterglow 
light curves were explored. In Paper H (Zeh, Klose, & Kann 2005) the optical light curve 
parameters were derived for a large (and basically complete) sample of afterglows observed 
by the end of 2004, up to the launch of the Swift mission. In this paper, we extend this 
systematic analysis to the spectral energy distribution (SED) of GRB afterglows in the 
optical/near-infrared bands in order to search for signals from extinction by dust in the 
GRB host galaxies. 

In § 2, we present the methods with which we analyzed the afterglow SEDs. Among the 
59 GRBs studied in Paper H, 30 had data quality sufficiently high to be included in a sample 
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for an investigation of the SEDs. We then further reduce this sample to a Golden Sample of 
19 GRB afterglows with parameters derived from specific dust model fits. In § 3 we present 
the results derived from our analysis and discuss our findings in the context of the standard 
fireball model. The fitting process allows us to derive the host galaxy extinction Ay along 
the line of sight and the intrinsic, extinction-corrected spectral slope (3 {F^ (x of the 
afterglow light in the optical/near-infrared bands. The Ay-f^ sample is then further analyzed 
statistically to derive conclusions about the environment of GRBs and the dust properties 
of high-redshift galaxies. 


2. DATA ANALYSIS 
2.1. The fitting Procedure 

In Paper II we presented the analysis of a complete list of optical/NIR afterglow light 
curves for which sufficient public data was available, up to GRB 041006. One parameter 
derived in these light curve fits is the normalization constant m^, the magnitude of the 
afterglow at a certain time, being one day after the GRB trigger in the case of a light curve 
fit with a simple power-law (iT(t) oc t““), or the jet break time, tf,, in case of a fit with a 
smoothly broken power-law with a pre-break decay slope oi and a post-break decay slope 
a 2 (Beuermann et al. 1999). All data used to fit these light curves is corrected for Galactic 
extinction with E{B — V) derived from the GOBE maps (Schlegel, Finkbeiner & Davis 
1998) employing the Milky Way extinction curve of Gardelli, Glayton & Mathis (1989). The 
contemporaneous afterglow brightness in different photometric bands is then transformed 
into a spectral energy distribution, using zero fluxes and median wavelengths taken from 
the literature. The SED for GRB 030329 is derived in an alternate way, see Appendix B. 
The fits use a Levenberg-Marquardt minimization algorithm. Unless stated otherwise, all 
errors in this paper are at the Icr confidence level. 

Initially, we assume no extinction and the SEDs are fit with a simple power-law, Fj, oc 
In the following, we label a slope derived from such a fit /5o. A steep /5o combined with 
a non-linear (curved) Fy is then indicative of extinction in the GRB host galaxy. 

In order to derive the visual extinction Ay in the GRB host galaxy along the line of 
sight, we transform the wavelengths of the SED into the host galaxy frame following (Fynbo 
et al. 2001b), using redshifts taken from the compilation of Friedman & Bloom (2005). In 
two cases redshift estimates are used. For GRB 980519, we use a redshift of z = 1.5, following 
Jaunsen et al. (2001) who state that .2 > 1.5 due to the lack of a supernova bump. For XRF 
030723, we varied the redshift and fit the SED for each redshift, and find the best results 
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in a redshift range oi z = 0.3 to 0.4. Thus, we adopt z = 0.35. Fynbo et al. (2004) find 
tentative evidence for a low redshift (although the host galaxy would be very faint), and 
Butler et al. (2005a) find a preferred 2 ; = 0.4 from the Amati relation (Amati et al. 2002). 
In both cases, the SED is fit well by a power-law with an additional small amount of source 
frame extinction. The influence of intergalactic Lyman absorption on the photometry was 
accounted for by excluding SED data points lying beyond 2.4 • 10^® Hz in the host frame 
from the fith The observed SED (corrected for Galactic extinction) is then modeled by the 
function 

F; = To • (1) 


with ^ 

I.U80 


( 2 ) 


Here, (3 is the intrinsic power-law slope of the SED, and Fq a normalization constant (we 
choose the unextinguished flux density at 5-10^^ Hz in the host galaxy frame). The function 
vi^host) = ^Ahost/^v is the extinction law assumed for the interstellar medium (ISM) of the 
GRB host galaxy. We call this extinction source frame extinction. It encompasses local 
extinction close to the site of the GRB and host galaxy extinction further away in case the 
afterglow passes through a significant part of the host galaxy along the line of sight. The 
extinction law for the Milky Way (MW), the Large Magellanic Gloud (LMG), and the Small 
Magellanic Gloud (SMG) was taken from Pei (1992). These three dust types differ strongly in 
the wavelength region we examine (host frame UV/optical/NIR), especially toward shorter 
wavelengths. The sequence MW, LMG, SMC is given by a decreasing strength of the 2175 
A bump, an increasing FUV extinction, and a decreasing reddening per H atom (Draine 
2000). The former implies the absence of very small carbonaceous grains (Weingartner & 
Draine 2001), while the latter is consistent with the trend observed for the metallicity in HII 
regions, with MW representing the highest metallicity and SMC the lowest. Thus, SMC dust 
traces the dust properties in low-metallicity environments. As the wavelengths of the SED 
are shifted from the UV into the optical/NIR, a distinction between the three dust types is 
easiest for GRBs at high redshift, whereas for nearby GRBs the three extinction curves are 
almost identical in the considered wavelength regime. 


Since the fit has three free parameters (Fq, [3 and Ay), the observed SED has to have at 
least four data points. The one exception is GRB 000131 (Appendix B), which we include 
because of its very high redshift {z = 4.5, Andersen et al. 2000). The fitting process is 
very sensitive to slight variations in the SED, especially if the data points have small 


^For GRB 030323 {z = 3.37), the V band also had to be excluded as the host galaxy is a DLA with a 
very wide Lyman a absorption line (Vreeswijk et al. 2004). 
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errors, so we did a careful check for each SED and removed outliers. Furthermore, we added 
a systematic error of 0.03 magnitudes to each data point of the SED (with the exception 
of those small number of cases where is based on only one data point, as these almost 
always have much larger errors anyway) to encompass uncertainties that may derive from 
the process of data reduction. This results in a strong decrease of for the fits and reduces 
the importance of this measure for discerning among different models for SED fits. 


2.2. Selection of a Golden Sample 

Our final sample (Table 1) comprises 30 afterglows out of the 59 bursts included in the 
sample of Paper II, with fits to the three dust models for each afterglow. Figure 1 shows all 
30 SEDs and the associated best fits for each of the three dust models. This figure contains 
all fit results that we obtain, including those with unphysical results, just in order to show 
the failure of certain extinction laws to reproduce the observational data. 

Table 1 reveals that many of the fits are unsatisfactory. First of all, some fits have very 
large error bars ATy and A/3. Since the SEDs are fit by allowing all three free parameters 
(inch Fo) to vary simultaneously, this creates a cross-sectional error, as the errors of [3 and 
Ty are correlated. The steepest spectral slope (/3 + A/3) results in the smallest extinction 
(Ty — Ady) and vice versa. Secondly, many fits, especially when using the MW dust model, 
find negative extinction. This would imply bright emission features, for which we have no 
evidence. Thirdly in some cases /3 < 0 is found, while the standard fireball model implies 
/3 > 0 in the timespans (hours to days) and wavelength bands that afterglows are typically 
observed in, and thus the timespan where we derive m^. These afterglow SEDs (e.g. GRB 
971214) are categorized by a relatively shallow unextinguished slope /3o ~ 0.5 to 1 but a 
strong spectral curvature. Correction of this curvature results in /3 < 0. One solution 
would be dust with an increased FUV extinction, but recent results on high redshift quasars 
(Maiolino et al. 2004) do not support this idea. We note that in these cases the data base is 
also sparse in photometric bands other than Rc- 

To derive a more homogenous sample, which we call the Golden Sample, we employ the 
following criteria: 

1. The Icr error in /3 (A/3) and the Icr error in Ay (ATy) should both be < 0.5. 

2. j4y -|- Aj4y ^ 0. 

3. We do not consider GRBs where all fits (MW, LMC, and SMC) find Ay < 0, even if 
the previous criterium is fulfilled. 


6 


4. > 0 (although we do not reject cases with (3 — A(3 < 0). 

5. A known redshift. 

As there is no case where LMC dust is clearly preferred compared to SMC and MW dust 
(Table 1), in the following we remove all LMC dust cases from further consideration. 

After applying these criteria, for the remaining 19 GRBs, in eight cases only one dust 
model remains. For GRB 000131 no distinction can be made, as it is not fit with extinction. 
The remaining ten GRBs give ambiguous results when it comes to distinguishing between 
the two dust models. With the exception of GRB 020813 and GRB 030328, where visual 
inspection of the fits shows that the inclusion of SMC dust reproduces the observed SED 
better than MW dust, no true distinction can be made at all. For these GRBs, we do not 
claim any dust model preference, but for the sake of consistency, we choose the dust model 
result which produces smaller errors in the derived parameters. This results in 19 P-Ay pairs 
which we designate our Golden Sample (Table 2). 


2.3. Highly extinguished GRBs outside our sample 

There are some GRBs for which high visual source frame extinction has been deduced in 
the literature but which are not included in Table 1. For GRB 980329, no redshift is known, 
although it is assumed to be high. Jaunsen et al. (2003) derive a Bayesian photometric 
redshift oi z = 3.6. The very steep observed spectral slope (uncorrected for extinction) has 
been noted by several authors. Yost et al. (2002) use their own data combined with data 
from the literature and derive an Ay ~ 1 for .2 = 3, corresponding to Nh 2 • lO^Ym ^ 
assuming a Galactic dust-to-gas ratio. This is concurrent with our result Ay = 1.03±0.65 for 
SMG dust and under the assumption of a redshift of ;2 = 3.6. However, this result also has 
/? < 0, implying that the true host extinction is lower than Ay — 1. Two afterglows which 
have redshifts but viable data in only two colors, GRB 990705 and GRB 000418, could not 
be fitted with the procedure outlined in §2.1. Masetti et al. (2000) find Pq = 2, in agreement 
with our result. While they do discuss reddening by intrinsic dust, their lack of a redshift 
measurement does not allow them to find a definite result. The steep initial decay suggests 
a wind environment with a cooling break blueward of the optical bands. Using this case and 
SMG dust, we derive Ay = 1.15. For GRB 000418, assuming a wind (ISM) environment with 
the cooling break blueward of the optical bands and the new derived light curve parameters 
from Paper II we fix P via the a-P relations and find Ay = 1.01 (Ay = 0.77). The latter is 
slightly smaller than in Klose et al. (2000) (who used an ISM environment), while Berger et 
al. (2001a) find Ay = 0.4 for LMG dust. Finally, there are indications of very high source 
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frame extinction from the NIR afterglows of GRB 030528 (Ran et al. 2004) and GRB 040827 
(De Luca et al. 2005), but we were not able to perform any fits. 


3. RESULTS AND DISCUSSION 
3.1. The prevalence of SMC-like dust in GRB host galaxies 

In the selection of our Golden Sample, there were eight bursts where application of our 
selection criteria yielded an unambiguous result. SMC dust is by far the most preferred 
model, with seven out of the eight cases (Table 2), only GRB 970508 has evidence for MW 
dust. This result is expected from studies of GRB host galaxies, which indicate that most 
hosts are low-metallicity galaxies, as they are Lyman a emitters, blue, and subluminous 
(e.g., Fynbo et al. 2003; Le Floc’h et al. 2003; Jakobsson et al. 2005b). Thus we deduce that 
the ISM in these galaxies is probably metal-poor. 

While the preference of SMC-like dust in GRB environments has already been found 
by many groups for several bursts, e.g., GRB 000301C (Jensen et al. 2001), GRB 000926 
(Fynbo et al. 2001b), GRB 010222 (Lee et al. 2001), GRB 020813 (Savaglio & Fall 2004), 
GRB 021004 (Holland et al. 2003), GRB 030226 (Klose et al. 2004), GRB 030429 (Jakobsson 
et al. 2004a), XRF 030723 (Fynbo et al. 2004), our study puts this conclusion on a statistical 
basis. Given the sparse knowledge we have so far on the dust properties in cosmologically 
remote galaxies we note however that this finding does not necessarily imply the requirement 
of low-metallicity environments for the creation of GRB progenitors, as it might be indicated 
by SMC-like dust. The fact that at least some afterglows SEDs are fitted best assuming MW 
dust (e.g., the spectral feature found in the spectrum of GRB 991216, Vreeswijk et al. 2005, 
which is mirrored in our SED result, and the strong preference of MW dust for GRB 030227, 
although the large errors make this case unsure) raises doubt on such a general requirement. 
Even though one can not rule out the possibility that in such cases the line of sight passes 
through foreground material with different grain properties and ISM metallicity. In any case, 
it is clear from this study that GRB afterglows can be used as a tool in order to explore 
the properties of cosmic dust in the cosmologically remote universe, in star-forming regions 
in particular and in galaxies in general. But high quality photometric data are essential for 
building larger samples. 
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3.2. The extinction within GRB host galaxies 

Among the 19 GRBs in our Golden Sample (Table 2), we find eleven with Ay — ^A-M > 0. 
We note that almost all of these cases are also those with the highest quality data. Basically, 
this is not surprising, as more SED data reduces errors. It is more speculative to state that, 
given high quality SED data, there will always be small amounts of source frame extinction 
detected in GRB afterglows. In eight of these cases, the detection is significant at the 2a 
level or higher, it is significant at more than 3cr in two cases and at 4cr in one case. 

The distribution of derived intrinsic extinctions is given in Eigure 2. A strong clustering 
toward low extinctions (Ay < 0.2) is evident. Zooming in with smaller bins shows that 
values between Ay = 0.1 and Ay = 0.2 are preferred, with a non-error-weighted mean value 
of Ay = 0.21 ± 0.04. No afterglow is extinguished by more than 0.8 magnitudes. As a 
comparison, the prime example of a dark burst is plotted, GRB 970828. For this burst, 
Djorgovski et al. (2001) derive Ay > 3.8 by interpolating between X-ray and radio data and 
comparing this with deep optical observations that found no afterglow Groot et al. (1998). 
The intervening space is free of examples, creating a potential Dark Burst Desert. On the 
other hand, there are two afterglows in Table 1 that have mean Ay values lying in this Dark 
Burst Desert but are not included in the Golden Sample due to the large errors of their 
fits (criterium (1) in § 2.2). For GRB 010921, MW dust is preferred, as SMG dust finds 
/? < 0. We derive Ay = 0.91 ± 0.82, in agreement with Price et al. (2003). While this is a 
positive detection, the error encompasses most other bursts of the Golden Sample. For GRB 
980703, no dust model preference can be derived, but even for the SMG result with smaller 
extinction, we find Ay = 1.32 ± 0.59. 

How real is the non-existence of highly extinguished afterglows? A first argument against 
the reality of the Dark Burst Desert is that our sample is biased toward afterglows with low 
source frame extinction, as we require viable photometric data in several wavebands and a 
known redshift. Intrinsic extinction complicates the gathering of photometry (starting with 
the discovery of the afterglow) and a spectroscopic redshift. In § 2.3, we mentioned several 
probably highly extinguished afterglows that are not part of our sample due to the sparsity 
of their photometric data. Highly extinguished afterglows are not observed well enough 
to determine conclusively that they are highly extinguished and not just intrinsically faint 
(Klose et al. 2003). 

Figure 3 shows the derived visual extinctions (Table 2) plotted as a function of redshift 2 ; 
(values taken from Friedman & Bloom 2005). With the exception of GRB 030429 the amount 
of extinction drops sharply toward high redshifts. Most likely, this can be explained by an 
observational bias, at least in the pre-Swift era where only very few afterglows were already 
localized within minutes after a burst. In particular, due to the wavelength dependence 
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of the dust opacity, the higher the redshift the more effective dust can dim an afterglow 
in the optical bands. Furthermore, the redshift measurement itself is biased toward bright 
and unextinguished afterglows, as many redshifts are determined from absorption spectra 
taken when the afterglow is much brighter than the underlying host galaxy. Finally, our 
data indicate that, on average, low-;^ GRBs have brighter afterglows (cf. Fignre 7), and thus 
have a greater chance of being detected through significant extinction. Unfortnnately, the 
chances of detecting a ’’Rosetta Stone Dark Burst” at low redshift that is highly obscured 
by dust but still bright enough to yield viable data and a redshift might be small if the GRB 
freqnency is coupled with the star-formation rate (e.g., Lloyd-Ronning et al. 2002; Firmani 
et al. 2004) and low metallicities (MacFadyen &: Woosley 1999). 

It is also visible that the the ability to discern between different dust models is coupled 
to the redshift, as the strongest deviations (2175 A bump, FUV extinction) appear in the 
UV region and the burst has to lie at a certain redshift to shift these features into the optical 
bands. 


3.3. Host extinction versus star formation rate 

Assuming that our database is not seriously affected by an observational bias (dark 
bursts). Figure 4 shows that the standard GRB afterglow is nearly unextinguished in its 
host galaxy. Since any long-lasting dust destruction by the intense fireball radiation shonld 
represent itself in a color evolution of the afterglow (Waxman & Draine 2000), which we do 
not find in our data, we conclude that on average there is not much evidence for dust in 
GRB environments along the line of sight. Furthermore, any dust destruction could only 
affect the dust in the immediate GRB environment, say, within 10 pc around the burster 
(Waxman &: Draine 2000; Draine 2000; Rhoads & Fruchter 2001; Draine & Hao 2002). Dnst 
at larger distances will still produce extinction. Additionally, the Swift satellite has made 
observations of early GRB afterglows routine, and a major result is the almost complete 
absence of reverse shock UV hashes (Roming et al. 2005). It seems that in most cases the 
absence of large amounts of dust is not dne to dust destruction but due to the fact that 
the standard GRB progenitor (seen pole-on along its rotational axis) is not enshrouded by 
dust, globally and locally. Thus, if in those cases where significant extinction is found the 
extinction is not just produced in the immediate GRB environment (local extinction) but 
prevails within the host galaxy (global extinction), we would expect that GRBs for which 
we find significant optical extinction are located in very dusty, submm-bright hosts. 

In order to explore this possibility, we have taken from the literatnre all snbmm data 
for GRB hosts with a flux density Fgso > 0 at the Icr level that are also in our main sample 
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of 30 afterglows (Berger et al. 2003a; Tanvir et al. 2004) and calculated the corresponding 
star formation rate (SFR) via the procedure developed by Yun & Carilli (2002) (cf. Berger 
et al. 2003a) to take the different redshifts into account. Since the SFR is a direct measure 
of the far-infrared luminosity of a galaxy (Kennicutt 1998), it traces its total amount of 
radiating dust. For GRB 980703, we conservatively used the result for SMC dust, which 
has the lowest extinction, and took the SFR from Berger et al. (2001b). We added GRB 
000418 {Ay = 1.01, § 2.3), arbitrarily assuming a 1 cr error of 0.5 mag, and transformed the 
Ay result for GRB 011211 (SMC dust) into an upper limit, as we find /3 < 0, making Ay 
unsure. The resulting relation of SFR (Mq per year) versus Ay is shown in Figure 4. A 
trend is visible, a rising SFR is coupled to a rising intrinsic visual extinction along the line of 
sight, with the exception of GRB 010222. A linear fit to the data (excepting GRB 010222^) 
gives SFR (Mq/yv) = 220 ± 130 -|- (500 ± 180) • Ay. We conclude that, on average, GRB 
afterglows that show significant extinction along the line of sight in their host are located 
in galaxies with a substantial star formation rate and, hence, a globally acting extinction 
by large amounts of dust. If this interpretation is correct, then the trend seen in the data 
indicates that the line-of-sight can pass through a significant extent of the host galaxy. Since 
the sample is still very small and the error bars on both the SFR and the extinction Ay are 
very large, it is clear that more and better data are required in order to verify this result. 


3.4. The intrinsic spectral slope and the power-law index of the electron 

distribution function 

Figure 5 displays the distribution of the intrinsic spectral slopes, (3, for the bursts from 
Table 2. The distribution is broad, ranging from 0.2 to 1.2, and features a peak around 
(3 = 0.7. We find a non-error-weighted mean value of (3 = 0.57 ± 0.05. This result of a broad 
(3 distribution of the intrinsic SED of the afterglows is probably robust, since (1) our light 
curve fits of the individual afterglows usually include many data points covering several days, 
so that we are not much sensitive to individual measurement errors and (2) these fits do not 
include the very early phase of an afterglow when its spectral properties might develop much 
faster than at later times. In fact, with the exception of the early afterglows of GRB 021004 


^The most remarkable outlier in the potential SFR-dy relation is GRB 010222. A Hubble Space Telescope 
image of the GRB 010222 host galaxy reveals the location of the GRB being offset by a small margin from 
the center of the galaxy (Pruchter et al. 2001; Galama et al. 2003). We find a low extinction value along the 
line of sight, whereas persistent submillimeter flux indicates that this is a dusty starburst galaxy (Frail et 
al. 2002). The discrepancy is resolved if the GRB happened toward the edge of the galaxy, but our line of 
sight places it in front of the galaxy. 
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and GRB 030329, we have never found clear evidence for color variations of the genuine 
afterglow light in our data. 

While it is not the goal of the present paper to find reliable explanations for the observed 
width of the distribution of spectral slopes [3 of the afterglows in the optical/NIR bands, we 
note that the mean of this distribution gives a value for the power-law index of the electron 
distribution {N{E)dE oc E~^dE) of p = 2.4, assuming a wind environment and the cooling 
frequency blueward of the optical bands (Zhang & Meszaros 2004, and references therin) 
and the peak of the histogram (Figure 5). This is in general agreement with theoretical 
predictions for ultra-relativistic shocks (Kirk et ah 2000; Achterberg et ah 2001). The max¬ 
imum values we find for [3 (Figure 5) could then be explained by those afterglows which had 
the cooling frequency redward of the optical bands during the entire time span when they 
were observed, provided that p is a universal number. On the other hand, a universality 
of p seems to be difficult to reconcile with those afterglows which have (3 < 0.5. Assuming 
the standard afterglow models (Zhang & Meszaros 2004), these afterglows require p < 2, 
provided that the basic model assumptions are indeed fulfilled in these cases. Note that in 
Paper II we also found afterglows that require p < 2 based on their light curve shapes alone. 

The fact that some afterglows require p < 2 has already been reported and in detail 
explored by others (cf. Panaitescu & Kumar 2000; 2001; Panaitescu 2005). Based on our 
SED fits three of 19 afterglows have (3 < 0.5 within their Icr error bar (GRB 970508, GRB 
991216 and GRB 030429) which implies p < 2, while a very shallow afterglow decay slope 
02 brings GRB 990123, GRB 991216, GRB 010222, GRB 030328 and GRB 041006 into this 
sample (Paper II). Obviously, only a minority of afterglows require p < 2. It is not obvious 
why these afterglows are specific in some sense. The afterglow of GRB 991216 is one of the 
brightest afterglow for redshifts z < 1. On the other hand, the afterglow of GRB 041006 is 
one of the less luminous afterglows (§ 3.6). Thus, p < 2 is not a question of luminosity. Note 
that both afterglows have only small intrinsic extinction (Ay ~ 0.1), so these low values for 
[3 are not just an artifact of the fitting process. For GRB 991216, the unextinguished slope is 
/5o = 0.54±0.03, for GRB 041006 it is /5o = 0.49±0.05. We also note that all these afterglow 
light curves are sampled fairly (GRB 990123) to very well (GRB 010222), implying that the 
finding of a flat spectral slope is not a question of data quality either. 


3.5. Dust-to-gas ratios in GRB host galaxies 

X-ray observations or the modeling of Lyman a absorption allow the determination of 
the hydrogen column density Nh along the line of sight to the GRB in its host galaxy after 
correcting for the column density in our Galaxy (Dickey & Lockman 1990). In Figure 6 we 
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present a sample of afterglows for which we have derived Ay and for which Nh values are 
reported in the literature (Table 3). This is an update of the plot first presented in Galama 
& Wijers (2001) and expanded in Stratta et ah (2004) The results reinforce our findings 
presented in § 3.1. With two exceptions, all points in the plot lie on the theoretical prediction 
for SMC dust or even above it. The two exceptions are GRB 021004 and 030329, where only 
upper limits for Nh are given in the literature, leading to dust-to-gas ratios higher than even 
for the Milky Way. On the other hand, both afterglow SEDs are best fit with SMC dust 
(although the preference is only weak for GRB 030329), implying that the three dust models 
we use are not applicable in all cases. Incidentally, GRB 021004 and GRB 030329 have 
the two best observed afterglow light curves with the most pronounced substructure (Paper 
II). The strongest outlier in Figure 3 is GRB 990123, for which we find a very low source 
frame extinction. The extremely bright UV flash of this burst (Akerlof et ah 1999) may 
have burned significant amounts of dust along the line of sight (cf. Waxman & Draine 2000; 
Galama & Wijers 2001; Rhoads & Fruchter 2001; Perna, Lazzati &: Fiore 2003), reducing 
the dust-to-gas ratio. No early multicolor data exists to probe the color variations that are 
expected in the early light curve. 


3.6. The luminosity distribution of the afterglows 

Knowledge of the intrinsic spectral slope of the afterglows allows us to determine their 
luminosity distribution (Appendix A). 

At first, in Figure 7, we show the Rc-band light curves of all 30 afterglows of the SED 
sample (Table 1) plotted with smoothed splines connecting the data points to guide the eye 
(for reasons of clarity, error bars on the photometry have been omitted). The data have been 
corrected for Galactic extinction and host contribution. For the light curve of GRB 030329, 
the supernova contribution has also been subtracted, using the data from Zeh, Klose, &: 
Hartmann (2005). A very large spread of magnitudes is seen, the range at one day after the 
trigger is 7.5 magnitudes between the afterglows of GRB 021211 and GRB 030329. Further 
bright afterglows are those of GRB 991208 and GRB 991216, other faint afterglows are GRB 
040924, GRB 030227 and GRB 971214. 

Figure 8 displays the afterglow light curves of the Golden Sample (Table 2) after applying 
the cosmological ^-correction (the second term of eq. A5) (Table 4) and after the time shifting 
to a common redshift of = 1 (the first term of eq. A5). In other words. Figure 8 is a 


^Note that we, unlike Galama & Wijers (2001) and Stratta et al. (2004), use a log-log plot to avoid 
crowding toward low Ay, leading to a linear depiction of the dust model curves. 
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measure of the absolute Rc band magnitudes of the afterglows up to a constant. Compared 
to Figure 7 the magnitude range has decreased, now being 5.7 magnitudes at one day after 
the burst. Nine of the 16 afterglows that have data one day after the burst lie in a range only 
two magnitudes wide, approximately clustered around the afterglow of GRB 030329. In other 
words, this afterglow is now seen to be quite typical. The two afterglows above this range 
are those of GRB 021004 and GRB 991208 (assuming for the latter burst an extrapolation of 
the decay with a ~ 2.5 to one day after the burst in the host rest frame. Paper II). The six 
afterglows beneath this range are GRB 030226, GRB 020405, GRB 030328, GRB 011121, 
GRB 041006 and GRB 040924. Of these six afterglows, only two (GRB 030226 and GRB 
030328) lie at z > 1. GRB 021004 is the most luminous afterglow at all times, although it 
is possible that the afterglow of GRB 991208 was brighter at earlier times when it was not 
yet discovered (Gastro-Tirado et al. 2001)^. 

Knowledge of the unextinguished light curves allows us to determine the luminosity 
distribution of our afterglow sample. Based on Fig. 8 we derive the Rc band magnitudes 
at one and four days (corresponding to half a day and two days after the burst in the host 
frame at z=l) and transform them into absolute magnitudes Mb. We do not extrapolate 
the afterglow light curves, except for GRB 030429 and GRB 020813, which both have their 
final data points close to 4 days and are post-break. Thus, not all fight curves are included 
(e.g., GRB 991208 is not yet detected at one day, GRB 030328 is not detected any more 
at four days if at z = 1). The sample then contains 16 GRBs at both one day and at four 
days, the results are given in Table 5. Thereby, we computed the luminosity errors in a 
conservative fashion. Given that the errors of (3 and Ty are correlated (the result Ty + ATy 
is coupled io (3 — A(3 and the other way around), we calculate three values of Mb for the pairs 
{13, Ty), (/? A/5, Ay — ATy), and (/5 — A/5, Ay + ATy). The pair (/5 A/5, Ay — ATy) 

results in a lower luminosity, and the pair (/5 — A/5, Ay + ATy) in a higher luminosity. The 
difference between these luminosities and the luminosity derived from the pair (/5, Ay) is a 
conservative upper limit for the error of Mb. In addition, we impose limits {3 — A/5 > 0 and 
Ay — ATy >0. If /5 — A/5 < 0 or Ay — ATy < 0, we set /5 = 0 or Ay = 0, respectively, 
when computing Mb. 

The resulting luminosity distributions are given in Figures 9 and 10. In addition, in 
order to search for a potential evolutionary effect or an observational bias, we distinguish 
GRBs with z < lA and ^ > 1.4, this being the median of the redshift distribution of our 
sample. A bimodal distribution of the afterglow luminosities is evident in Figures 11 and 
12. Afterglows sA, z < 1.4 tend to be less luminous, which might be explainable by an 


"‘We note that both Nardini et al. (2005) and Liang & Zhang (2005) have reached similar conclusions 
using host galaxy extinction values derived from the literature. 
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observational bias: the chances to detect intrinsically faint afterglows are higher for a lower 

z. 


At one day after the GRB at z = 1, the complete distribution has a two peaks, the 
mean lying at Mb = —23.2 ± 0.4, with the most luminous afterglow, GRB 021004, having 
Mb = —25.58lo]o45 ^^cl the faintest, GRB 040924, having Mb = — 19.84]'^Q;g5. The means 
for the distributions of the z < 1. 4 and .2 > 1.4 GRBs are Mb = —22.4 ± 0.6 and Mb = 
—24.1 ±0.5, respectively. The difference between the mean values is 1.7T0.7 mag, indicating 
a bimodality. 

At four days after the burst at z = 1, a single broad peak remains, for the complete 
distribution, the mean value is Mb = —21.4±0.4. Now, the brightest afterglow (GRB 021004) 
has Mb = —24.07lo!o4 ^^e faintest afterglow (GRB 041006) has Mb = —18.89lo!65- 

Once again, the bimodal distribution for nearby {z < 1.4) and distant {z > 1.4) GRB 
afterglows is evident. The means for the distributions of the z < 1.4 and z > 1.4 GRBs are 
Mb = —21.2±0.5 and Mb = —21.8±0.6, respectively. The difference between the mean 
values is 0.6 ± 0.8 mag, the significance of the bimodality has been strongly reduced. 

Gendre & Boer (2005) analysed the X-ray afterglow light curves of GRBs and shifted 
them to 2 ; = 1 in an analog process, finding a bimodal flux distribution. There are eight 
GRBs in their sample (GRB 970508, GRB 990123, GRB 991216, GRB 000926, GRB 010222, 
GRB 011121, GRB 030226, and GRB 030329) that are also among the afterglows we have 
plotted in Figure 8. GRB 990123, GRB 991216, GRB 000926, and GRB 010222 are in 
group I of Gendre & Boer (2005), while GRB 970508, GRB 011121, GRB 030226, and GRB 
030329 are in their group 11. We note that iox z = 1, t = 1 day lies at the end of the initial 
plateau phase of the light curve of GRB 970508. The afterglow brightens by 1.5 magnitudes 
shortly afterward. A comparison reveals that while the two groups mix in the optical (the 
afterglow of GRB 030329 being more luminous than those of GRB 990123 and GRB 010222 
at one day after the burst at .2 = 1), the mean absolute magnitude of group I afterglows is 
Mb = —24.2 ± 0.4, while it is Mb = —22.7 ± 0.5 for group II afterglows, with the difference 
being 1.5 ± 0.7 mag. While this finding is intriguing, the sample is too small to draw any 
conclusions. 

Knowledge of the extinction corrected afterglow magnitudes allows us to compare the 
luminosities of the afterglows at a common point in the evolution of the jet, the jet break 
time. Ignoring any fine structure in the light curves, for all bursts in our Golden Sample 
with a jet break (Paper II) we computed the apparent Rc-band afterglow magnitudes at the 
time of the jet break from the Beuermann equation (Beuermann et al. 1999). We included 
the fit around the jet break of GRB 030329 but excluded GRB 021004, as the reality of the 
late break we found in Paper II is unclear due to the many rebrightening episodes. These 
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magnitudes were then converted to luminosities and normalized to the luminosity of the 
afterglow of GRB 990123 at the time of its jet break. 

We use these luminosity ratios to search for correlations between afterglow luminosity 
and parameters of the prompt 7 emission of the GRBs. We take the isotropic energy release 
Eiso and the source frame peak energy Ep from Friedman &: Bloom (2005), and the collimation 
corrected energies E^ from Paper 11. We do not find any correlations between the afterglow 
luminosity during the break time and these three parameters, with the absolute value of the 
correlation coefficient being smaller than 0.25 in all cases. 


4. SUMMARY & CONCLUSIONS 

We have presented a sample of 30 GRB afterglow spectral energy distributions in the 
optical/NIR bands which have been modeled with various dust extinction curves (Milky Way, 
Large Magellanic Gloud and Small Magellanic Gloud) to derive the source frame extinction, 
4v, intrinsic to the host galaxies and the spectral slope, (3, of the afterglows unaffected by 
any dust extinction. As all afterglows have been analyzed in a systematic way, the results are 
fully comparable, making this sample unique in terms of both size and consistency. For the 
further statistical study, we selected 19 afterglows, our Golden Sample, which have physically 
reasonable results and small error bars. 

The preferred dust models we find (§ 3.1) as well as the deduced source frame dust- 
to-gas ratios (§ 3.5) based on the inclusion of data taken from the literature, both indicate 
that the majority of GRBs we have investigated, covering the redshift range from 0.1 to 
4.5, occur in low-metallicity environments. The distribution that we have derived from 
these data (§ 3.2) highlights a sparsity of strongly extinguished afterglows, creating a Dark 
Burst Desert, even though it is unclear if the preference of low extinctions is more than an 
observational and sample selection bias. Our finding that most afterglows suffer from only 
low extinction in their hosts could indicate that afterglows are usually not obscured by dust 
close to the burster. One would then expect that the most extinguished afterglows are in 
fact located in globally dusty hosts. Indeed, we find weak evidence for a correlation between 
the submm flux of GRB host galaxies and the source frame extinction Ay- Although the 
statistical significance is low due to the small sample size and the large errors, this finding 
calls for a more thorough investigation. 

Knowledge of (3 and Ay allowed us then to correct the afterglow light curves for intrinsic 
extinction and to derive the true luminosity distribution of our afterglow sample at chosen 
times in the host galaxy frame (§ 3.6). We find that, on average, low -.2 afterglows are less 
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luminous than high-^ afterglows. The most likely explanation we have at hand for this 
finding is an observational bias against intrinsically faint afterglows at high redshifts. A 
bimodal distribution found by Gendre & Boer (2005) in similarly corrected X-ray afterglows 
is not clearly seen in the optical although, on average, GRBs with fainter X-ray afterglows 
also have fainter optical afterglows. Unfortunately, the available sample size is still too small 
to reach definite conclusions. A search for correlations between prompt emission parameters 
and the luminosity of the optical afterglows at jet break time has come up empty. 

Since our sample is exclusively composed of GRBs from the pve-Swift era, a similar study 
in a few years time on Swift-discovered GRB afterglows will shed light on the Dark Burst 
Desert and the true afterglow luminosity distribution by removing observational bias factors 
via rapid and highly precise GRB localizations. Already, Swift has lead to the discovery of 
very faint afterglows (e.g. GRB 050126, GRB 050607, Berger et al. 2005a; Rhoads 2005) 
including what may be the ’’darkest” burst ever, GRB 050412 (Kosugi et al. 2005; Jakobsson 
et al. 2005a). The recent discovery of the first afterglows of short GRBs (Hjorth et al. 2005; 
Fox et al. 2005; Covino et al. 2005; Berger et al. 2005b; Soderberg & Berger 2005) opens the 
possibility of finally making a comparison of the environment of the two different classes of 
GRBs. 
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A. Shifting the afterglows to a common redshift 

In the following, we consider the proper afterglow, i.e., cleaned from any underlying host 
component as well as corrected for Galactic extinction and extinction in the host galaxy. In 
cases where the late-time light curve is dominated by a supernova component, the data 
was removed due to a lack of knowledge concerning the supernova SED. We start with the 
expression for the flux density per unit frequency, F{i',t), of a time-dependent source that 
is shifted to a redshift Zi, which was originally observed at a redshift Zq (Dado, Dar & de 
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Riijula 2002): 
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Here, t and u are measured in the observers frame, is the luminosity distance. The 
exponential function corrects for a postulated extinction in the host galaxy at redshift zq, 
where r is the optical depth in the host at the observed frequency. 

When shifting an afterglow to a redshift zi, its apparent magnitude changes by an 
amount Am = m{z\) — m{zo), with 
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where <S'(A) is the wavelength-dependent hlter response function for the given photometric 
band. For an unabsorbed afterglow 


F{zo-, iz, t) = const. f{t) p ^ . (A3) 

We assume (3 = const. Since the light curve shape can be affected by a jet break and by 
rebrightening episodes (which are achromatic according to our present data base), we do not 
specify f{t). Using eq. (A3), it follows that in eq. (Al) 
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with i = t Any extinction correction is included in F(zi; p,t) in eq. (Al). 


(A4) 


After transforming flux density per unit frequency into flux density per unit wavelength 
and inserting eqs. (Al, A3, A4) into eq. (A2), we obtain 
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In a logarithmic plot (apparent magnitude vs. log time) the first term in eq. (A5) represents 
a shift of the observed light curve in time. The second term in eq. (A5) shifts the light curve 
along the magnitude axis; the light curve shape is not affected. 
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In our calculations we assumed a flat universe with matter density VLm = 0.27, cos¬ 
mological constant Qa = 0.73, and Hubble Constant Hq = 71 km s“^ Mpc“^ (Spergel 
et al. 2003). For S{X) we used the Alter functions of the VLT/FORSl Bessel filters (see 
http: // WWW. eso .org/instruments / fors 1 / filters .html). 


B. Comparison with the literature 

We have searched the literature for publications that also derive the source frame ex¬ 
tinction toward the gamma-ray-bursts listed in Table 1. In most cases, these results are in 
good agreement with our findings. 

GRB 970508 To our knowledge, Reichart (1998) was the first who discussed the ex¬ 
tinction in a GRB host galaxy derived via a fit of the observed SED in the optical bands. 
Analysing collected data from the literature, he finds for GRB 970508 Ay = 0.24lo!o8 (albeit 
for a slightly higher z, having left the redshift as a free parameter), concurrent with our 
result Ay = 0.38 ±0.11 for the preferred MW dust extinction model. 

GRB 971214 Halpern et al. (1998) find a strong spectral curvature in their data and 
infer a large E{B — R) 0.4 for an assumed redshift of .2 = 2 . As a fixed value for the 
spectral slope, they took the V — I color of the afterglow of GRB 970508. This strong spectral 
curvature is mirrored in our SED fit which finds a large extinction and /? < 0. Eits with a 
fixed (3 derived from the a-(3 relations (Zhang & Meszaros 2004) are strongly rejected. 

GRB 980519 We have found no references to intrinsic extinction in the literature. This 
makes our preferred SMC dust result of Ay = 0.22 ± 0.19, derived under the assumption of 
a redshift of z = 1.5, a novel result. 

GRB 980703 This burst has the highest extinction in our sample (Table 1). This result 
is not unprecedented, however. While finding different values, several authors reported on a 
high extinction at different times after the burst: Castro-Tirado et al. (1999) find Ay = 2.2 
0.9 days post-burst, Vreeswijk et al. (1999) find Ay = 1.5 ± 0.11 at 1.2 days, and Bloom 
et al. (1998) derive Ay = 0.9 ± 0.2 at 5.3 days. Our achromatic fit gives Ay = 1.93 ± 0.91 
for the MW dust model, close to the result from Castro-Tirado et al. (1999), while SMC 
dust gives Ay = 1.32 ± 0.59. The light curves are not sampled well enough to check if the 
discrepancies between the various authors are the result of chromatic changes of the SED, a 
bumpy structure, or a true decrease of line-of-sight extinction with time. 

GRB 990123 Galama et al. (1999) find a negative extinction for this burst and therefore 
fix it to zero, consistent with the low value we derive: Ay = 0.04 ± 0.05 for the preferred 


19 


SMC dust model. On the other hand, Savaglio et al. (2003) derive a very high extinction of 
Ay = 1.1 from metal column abundances. Given the high quality of the optical photometry, 
we can not reproduce such a high extinction. 

GRB 990510 Both, Stanek et al. (1999) and Beuermann et al. (1999) find a slight 
curvature in the SED. Removing the R-band data point, Stanek et al. (1999) find [3 = 
0.46 ± 0.08, comparable to our result, {3 = 0.30 ± 0.69 for SMC dust. Neither publication 
fits the SED with a dust model to derive the host extinction. 

GRB 991208 Almost all data on this afterglow (in BVRcIc) is from Castro-Tirado 
et al. (2001). They do not discuss intrinsic extinction in the host galaxy. Together with 
additional A'-band data from Bloom et al. (1999) and Chary et al. (2002), we find a high 
extinction of Ay = 0.80 ± 0.29 for MW dust, the weakly preferred model, and an excellent 
fit (x^/dof=0.20). 

GRB 991216 After correcting for significant Galactic foreground extinction, both Gar- 
navich et al. (2000) and Halpern et al. (2000) find a SED well approximated by a power-law, 
concurrent with no dust reddening, although these authors do not discuss intrinsic redden¬ 
ing. With a larger database and more colors, we find Ay = 0.13 ± 0.08 for the weakly 
preferred MW dust, still a low value. This is consistent with the value Ay = 0.16 ±0.02 that 
Vreeswijk et al. (2005) find from a low-resolution spectrum. They also detect a broad ab¬ 
sorption feature centered at 2360 A, which can be interpreted as a redshifted 2715 A bump, 
giving credence to our choice of MW dust as the preferred model. 

GRB 000131 Eor an SMG dust model, Andersen et al. (2000) find Ay = 0.155 ± 0.045. 
We are not able to perform a free fit as the SED has only three colors [Re and V are affected 
by Lyman dampening), but fitting a straight power-law shows the SED is consistent with 
no extinction (i.e., no SED curvature is seen). 

GRB 00030IG Jensen et al. (2001) find very low host extinction (Ay = 0.09 ± 0.04) 
for this burst using an SMC dust curve, and no solution for MW and LMC dust. Rhoads 
& Fruchter (2001) derive a best fit result of Ay = 0.09 for SMC dust, find a slightly worse 
fit for LMC dust and a much worse fit for MW dust. We also find negative extinction for 
MW dust, and also low values for LMC and SMC dust (Ay = 0.12 ± 0.06 for SMC dust), 
completely in concordance with the results of the other authors. 

GRB 000911 Comparing their host galaxy SED with synthetic spectra of extinguished 
galaxies, Masetti et al. (2005) find Ay = 0.32 for SMC dust, consistent with our result of 
Ay = 0.27 ± 0.32. The LMC result (Ay = 0.39 compared to Ay = 0.27 ± 0.30) is also in 
agreement. 
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GRB 000926 For SMC dust, we find Ay = 0.15 ± 0.07, in good agreement with Fynbo 
et al. (2001b) who find = 0.18 ± 0.06 for SMC dust. Several other analyses (Price et ah 
2001; Sagar et al. 2001) disagree more with our results, but Fynbo et al. (2001b) employ a 
UBVRcIcJHK spectral energy distribution as we do, while the other studies have shorter 
baselines. Based on metal column abundances from high-resolution spectra, Savaglio et al. 
(2003) find a large Ay = 0.9, which we can not confirm via the broad band photometry. 

GRB 010222 Working with five-color SDSS photometry, Lee et al. (2001) find low 
extinction (j4v < 0.057 for [3 = 0.75 fixed) for this burst with an SMC dust curve, which 
is less than what we derive {Ay = 0.14 ± 0.08 for (3 = 0.76 and SMC dust). Our result is 
consistent with that of Galama et al. (2003), who find Ay = 0.1 ±0.02. Savaglio et al. (2003) 
find a higher extinction of = 0.7 based on metal column abundances. 

GRB 010921 Not finding a supernova bump in this nearby GRB, Price et al. (2003) 
invoke Ay 1. While Zeh, Klose, &: Hartmann (2004) do find a weak SN bump and we 
derive the SED from this fit, we also find significant extinction of Ay = 0.91 ± 0.82 for MW 
dust, the preferred model, thus validating the analysis in Price et al. (2003). 

GRB 011121 Price et al. ( 2002 ) deduce a Galactic visual extinction of Ay = 1.16 ±0.25 
along the line of sight by fitting the uncorrected data points with a MW extinction law at 
z = 0, which translates into a Galactic E{B — V) = 0.37 ± 0.08, assuming Ry =3.1. They 
could not constrain the host extinction, however. Garnavich et al. (2003) find a Galactic 
E{B — V) = 0.43 ± 0.07 and a spectral slope of /5 = 0.66 ± 0.13 between about 0.4 and 
1.5 days after the burst; a possible extinction in the host galaxy is not considered. Greiner 
et al. (2003) adopt E{B — V) = 0.46 based on the COBE data (Schlegel, Einkbeiner & 
Davis 1998), find no evidence for host extinction, and (3 = 0.62 ± 0.05 at 2.5 days after the 
burst. They note, however, that the HI maps of Dickey & Lockman (1990) favor a lower 
Galactic extinction of Ay = 0.9 mag, i.e. E{B — V) = 0.29 (assuming Ry = 3.1). We 
refitted the combined data set taken from Greiner et al. (2003), adding carefully selected 
(early and HST data) data from Garnavich et al. (2003) and Price et al. ( 2002 ). In addition, 
we have now used E{B — V) = 0.29 deduced via the HI data (Dickey & Lockman 1990), 
in order to overcome the potential mixture of an unknown source frame extinction and an 
uncertain Galactic extinction along the line of sight. We omitted the i7-band data (Greiner 
et al. 2003), however, since they worsen the fit. In doing so, we now derive a moderate host 
extinction of Ay 0.4 mag (SMC dust). The MW dust model is not preferred, although 
the low redshift of the GRB makes a differentiation among the three dust models hard. 

GRB 011211 Both Holland et al. (2002) and Jakobsson et al. (2003) find very low 
extinction for this burst. Jakobsson et al. (2003) rule out MW dust and find that SMG dust 
gives the best fit with Ay = 0.08 ±0.08. Our result is that MW dust is ruled out, LMG dust 
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gives a low Ay with a very bad fit, and SMC dust gives a moderate Ay = 0.25 ± 0.06 but 
with /I = 0 ± 0.15. Thus the results are comparable in principle. 

GRB 020124 Hjorth et al. (2003b) find Ay < 0.2 by fixing (3 > 0.5 for SMC dust, 
deeming (3 < 0.5 unrealistic. Their SED is consistent with ours, showing a strong curvature. 
They additionally employ a synthetic S-band point derived from extrapolating the spectrum 
and find [3 = 0.31 ±0.43 from a free fit of SMC dust, still higher than what we derive 
(excluding the B band, as it is affected by Lyman dampening), although the results are 
identical within the large errors. Berger et al. (2002) find Ay = 0 to 0.9. 

GRB 020405 Working with almost the same data set as we do, and including X-ray 
observations, Masetti et al. (2003) find negligible (though unspecified) dust extinction in this 
afterglow. On the other hand, Stratta et al. (2005), comparing the X-ray and the optical 
brightness, invoke a gray extinction curve with a very high Ay = 2.3. We find moderate 
extinction Ay ss 0.25 for all three dust models. 

GRB 020813 Savaglio & Fall (2004) performed a very stringent analysis of an early 
high-resolution spectrum of this afterglow and, via analysis of metal column densities, derive 
Ay < 0.08, Ay < 0.19 and Ay < 0.18 for MW, LMC, and SMC dust, respectively. Our 
analysis gives a very good match, with definite results: Ay = 0.01 ± 0.08, Ay = 0.19 ± 0.12 
and Ay = 0.12 ± 0.07, respectively. Furthermore, we also come to the conclusion that SMC 
dust is the preferred model. To make their unextincted spectral slope more compatible with 
standard models, Li et al. (2003) invoke E{B — V) = 0.05 (which converts to Ay = 0.15). 
Covino et al. (2003) deduce Ay — 0.12 ± 0.04, but for the Galactic extinction curve of 
Cardelli, Clayton & Mathis (1989). 

GRB 021004 For an SMC dust curve, Holland et al. (2003) derive Ay = 0.26 ± 0.04, 
coupled with (3 = 0 .39±0.12. This work finds a lower Ay = 0.14±0.05 coupled with a higher 
(3 = 0.67 ± 0.14. This is for the free fit, a wind model with the cooling break blueward of 
the optical gives [3 = 0.38 fixed and Ay = 0.20 ± 0.02, which is comparable. Fitting an SED 
derived on October 11.67, Fynbo et al. (2005) find exactly this result (ffy = 0.20 ±0.02) for 
SMC dust, associated with [3 = 0.42 ± 0.06. They find higher extinction for LMC dust and 
a worse fit, and negative extinction for MW dust, just as we do. 

GRB 021211 Fox et al. (2003) derive extinction values for this GRB by fixing (3. As our 
light curve fitting finds a very blue afterglow with f3o = 0.10 ± 0.09 (i.e., almost fiat), these 
values can not be compared adequately. Fox et al. (2003) find up to Ay = 0.64, while we 
only derive small negative extinction for all dust models. Since these authors have derived 
their spectral slope ~ 1 from an early B — Kg color, this may imply color evolution (cf. 
Nysewander et al. 2005). Holland et al. (2004) derive an unextincted slope of (3o = 0.69±0.14 
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at 0.87 days after the burst in the observer frame and an upper limit on the source frame 
extinction of Ay < 0.08. 

GRB 030226 Klose et ah (2004) find negligible extinction for this afterglow based on 
their multicolor data set, in contrast to Pandey et ah (2004) who, on theoretical grounds, 
argue for a fixed f3 = 0.55 and thus derive small extinction from their (Oq = 1. We utilize 
both data sets, and find a resnlt concurrent with Klose et ah (2004), MW and LMC dust are 
ruled out, and the reddening using SMC dust (j4v = 0.06 ± 0.06) is negligible within errors. 

GRB 030227 Castro-Tirado et ah (2003) were not able to conclusively fit their SED. 
Their (3q = 1.25±0.14 is in perfect agreement with the result we find, (Oq = 1.24±0.13. While 
the data set is sparse, it is, with the exception of one data point, from Castro-Tirado et ah 
(2003) exclusively. The steep /3o supports dust reddening, and the unextinguished (3 = 0.78, 
a reasonable value, further supports our findings of Ay = 0.38 ± 1.81, which could explain 
the faintness of this not too distant (z = 1.39, Watson et ah 2003) afterglow. The sparse 
data points of the SED and their large errors lead to very large errors in (3 and Ty, however. 

GRB 030323 The SED of this burst is unusual. The B band can not be included 
due to Lyman dampening, and the spectrum (a strong DLA) shows that the V band is 
also unreliable. With the exception of the K band, the SED shows a strong curvature 
without being very steep. As there is no reason to remove the K band point (based on 
three photometric points, and not two each as in J and i7), we tried to fit the SED with it 
(removing it results in very good fits that have /5 < 0). As the data points have small errors, 
the result is a very bad fit {xl ~ 6 for two degrees of freedom). Vreeswijk et al. (2004) do 
not fit the light curves, but analyse the SED at different post-bnrst times from observations 
at near-identical epochs. They fix /? = 0.28, derived from the a-(3 relations. Using this 
value, they find Ay = 0.5, Ay = 0.25 and Ay = 0.16 for MW, LMC and SMC dust curves, 
respectively. Eixing (3 to this value and fitting our SED yields Ay = 0.70, Ay = 0.41 and 
Ay = 0.26 for MW, LMC and SMC dust curves, respectively. While these values are higher 
than those in Vreeswijk et al. (2004), they are comparable if conservative errors are assumed. 
However, the A'-band data is a strong outlier in these fits. 

GRB 030328 After Paper II went into press, Maiorano et al. (2006) presented a large 
amount of UBVRcic data on this burst that was not available to us before. We refitted 
the optical light curves and derive the following parameters which supersede those presented 
in Table 1 of Paper II: x^/dof = 0.72 (1.34 in Paper II), dof = 81 (18 in Paper II), rric = 
19.54 ± 0.31 mag (20.61 ± 0.23 in Paper II), ai = 0.61 ± 0.13 (0.87 ± 0.04 in Paper II), 
= 1.41 ± 0.12 (1.54 ±0.11 in Paper II), 4 = 0.29 ± 0.06 days after the burst (0.60 ± 0.10 
in Paper II), n = 3.19±2.74 (fixed to 10 in Paper II). The host galaxy magnitude remains 
unchanged. (See Paper II for the definitions of these values.) These changes do not influence 
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the conclusions of Paper 11. We note that the value for the break smoothness parameter n 
is fully in agreement with the possible correlation between «i and n (Figure 8 of Paper II), 
and is almost identical to the value pair derived for GRB 010222. 

From the Fe II column density derived from the optical spectrum, Maiorano et ah (2006) 
estimate Ay <0.1 mag, in agreement with our result Ay = 0.05 ± 0.15 mag. 

GRB 030329 The SED of this burst was derived in an alternate way compared to all 
other SEDs. The rebrightening episodes starting at one day are found to be achromatic over 
several days. The high data quality makes it possible to shift the light curves to a common 
magnitude, the amount of shift is used to construct the SED, assuming conservative errors. 
This method is independent of the light curve fit, which is very complicated in the case 
of GRB 030329 (cf. Paper II). Bloom et al. (2004) find Ay = 0.94 ± 0.24, coupled with 
[3 = 0.11 ±0.22, for MW dust. Our result is: Ay = 0.54±0.22, coupled with (3 = 0.30±0.22, 
which is comparable within errors. Eor f3 = 0.5 fixed, Bloom et al. (2004) derive Ay < 0.3, 
comparable to Ay = 0.34 ± 0.04 which we derive when fixing (3 at this value. 

GRB 030429 Jakobsson et al. (2004a) rule out MW dust (as do we) and find Ay = 
0.30 ± 0.06 for SMG dust at 0.548 days after the burst, slightly lower than our value of 
Ay = 0.40 ± 0.10. They also find higher extinction for LMG dust, but the fit is much worse, 
contrary in part to our result. We find much lower source frame extinction for LMG dust, 
but here also the quality of the fit is bad. 

XRF 030723 Fynbo et al. (2004) give Ay < 0.5 for SMC dust and z = 0.3, completely 
concurrent with our result of Ay = 0.28 ± 0.24 for SMC dust and under the assumption of 
a redshift of ^ = 0.35. 

GRB 040924 Silvey et al. (2004) use their own data and data from the GCN archives 
to derive Po = 0.61 ± 0.08, slightly lower than the value we find (/3o = 0.80 ± 0.03). They do 
not fit the SED with a dust model. Soderberg et al. (2006) derive (3o ~ 0.7 and argue that 
for the standard blastwave model, (3 > 0.5 and thus Ay < 0.16 mag. This is in agreement 
with our result Ay = 0.16 ± 0.44 mag. 

GRB O 4 IOO 6 Soderberg et al. (2006) find /5o ~ 0.5 and, to keep p >2, assume negligible 
host galaxy extinction. This is in agreement with our result Ay = 0.11 ± 0.23, but we note 
that the shape of the SED is indicative of a slight amount of dust. The derived (3 = 0.36±0.27 
and the shallow post-break light curve decay slope 0:2 = 1.30 ± 0.02 (Paper II) are both 
indicative of a hard electron spectrum with p < 2. 
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Table 1. Results of the SED fitting 


GRB 


Filters^ 

^dof 

MW Dust 

n 

Av 

^dof 

LMC Dust 

0 


^dof 

SMC Dust 

0 


970508 


BVRcIcKs 

2.73 

0.32 ±0.15 

0.38 ±0.11 

2.00 

0.11 ±0.20 

0.55 ±0.15 

1.14 

0.00 ±0.22 

0.61 ±0.15 

971214 


RcIcJK 

0.32 

2.26 ±0.81 

-1.44 ±0.72 

0.14 

-1.26 ±0.95 

1.06 ±0.52 

0.09 

-0.50 ±0.58 

0.44 ± 0.21 

980519^^ 

1 

UBVRcIc 

0.61 

1.11 ±0.11 

-0.03 ±0.05 

0.73 

1.10 ±0.26 

-0.02 ±0.12 

0.01 

0.44 ± 0.54 

0.22 ±0.19 

980703 


BVRcIcJHK 

2.28 

0.57 ±0.92 

1.93 ±0.91 

2.10 

0.57 ±0.85 

1.85 ±0.81 

2.33 

1.05 ±0.66 

1.32 ±0.59 

990123 

2 

UBVRcIcHK 

0.21 

0.62 ± 0.09 

-0.06 ±0.07 

0.35 

0.50 ±0.20 

0.04 ±0.13 

0.20 

0.46 ±0.12 

0.04 ± 0.05 

990510 

3 

BVRcIc 

0.03 

0.89 ±0.11 

-0.05 ± 0.04 

0.01 

1.12 ±0.31 

-0.15 ±0.15 

0.58 

0.30 ±0.69 

0.18 ±0.24 

991208 


BVRcIcK 

0.20 

0.23 ±0.37 

0.80 ±0.29 

0.17 

0.07 ±0.43 

0.93 ±0.34 

0.26 

0.19 ±0.40 

0.76 ± 0.28 

991216 

4 

BVRcIcJHK 

0.12 

0.38 ±0.11 

0.13 ±0.08 

0.16 

0.32 ±0.15 

0.18 ±0.11 

0.29 

0.30 ±0.17 

0.18 ±0.13 

000131^ 


IcHKs 

<0.01 

0.66 ± 0.34 

0 

<0.01 

0.66 ± 0.34 

0 

<0.01 

0.66 ± 0.34 

0 

000301G 

5 

BVRcIcJK 

2.48 

0.88 ± 0.07 

-0.03 ±0.05 

2.03 

0.59 ±0.19 

0.16 ±0.11 

1.14 

0.59 ±0.12 

0.12 ±0.06 

000911 


BVRcIcJKs 

0.26 

0.75 ± 0.26 

0.20 ±0.22 

0.28 

0.67 ±0.36 

0.27 ±0.30 

0.33 

0.65 ± 0.40 

0.27 ±0.32 

000926 

6 

BVRcIcJHK 

1.02 

1.43 ±0.07 

0.07 ±0.05 

1.67 

1.29 ±0.20 

0.04 ±0.12 

0.37 

1.01 ±0.16 

0.15 ±0.07 

010222 

7 

UBVRcIcJK 

1.35 

1.18 ±0.07 

-0.03 ± 0.04 

1.36 

1.02 ±0.20 

0.07 ±0.10 

0.58 

0.76 ± 0.22 

0.14 ±0.08 

010921 


UBVRcr*i* 

0.07 

0.81 ± 1.21 

0.91 ±0.82 

0.06 

0.03 ± 1.88 

1.44 ± 1.29 

0.03 

-1.01 ±2.71 

1.91 ± 1.66 

011121 

8 

UBVRcIcJK 

0.55 

0.55 ±0.15 

0.49 ±0.18 

0.54 

0.55 ±0.15 

0.47 ±0.17 

0.49 

0.61 ±0.13 

0.39 ±0.14 

011211 

9 

BVRcIcJK 

2.03 

0.88 ± 0.09 

-0.24 ±0.06 

7.12 

0.38 ± 0.24 

0.11 ±0.14 

1.89 

0.00 ±0.15 

0.25 ± 0.06 

020124 

10 

RcIcJKs 

0.14 

1.29 ±0.65 

-0.44 ±0.59 

0.15 

-0.49 ± 1.77 

1.34 ±0.98 

0.01 

0.11 ±0.85 

0.28 ±0.33 

020405 

11 

BVRcIcJHKs 

1.82 

0.96 ± 0.20 

0.14 ±0.18 

1.75 

0.91 ±0.22 

0.19 ±0.20 

1.75 

0.94 ±0.19 

0.15 ±0.16 

020813 

12 

UBVRcIcJHK 

1.62 

1.03 ±0.11 

0.01 ± 0.08 

1.10 

0.74 ± 0.20 

0.19 ±0.12 

0.95 

0.81 ±0.14 

0.12 ±0.07 

021004 


BVRcIcJHKs 

1.62 

1.31 ±0.10 

-0.20 ±0.08 

1.57 

0.56 ±0.22 

0.27 ±0.11 

0.69 

0.67 ±0.14 

0.14 ±0.05 

021211 


BVRcIcJHKs 

1.16 

0.15 ±0.23 

-0.04 ±0.15 

1.13 

0.22 ± 0.32 

-0.09 ± 0.23 

1.05 

0.38 ± 0.40 

-0.19 ±0.27 

030226 

13 

BVRcIcJHK 

0.15 

0.77 ±0.06 

-0.06 ±0.04 

0.57 

0.78 ±0.18 

-0.05 ±0.11 

0.32 

0.57 ±0.12 

0.06 ± 0.06 

030227 


BRcHK 

1.41 

0.78 ±2.17 

0.38 ± 1.81 

0.74 

2.23 ± 1.15 

-0.76 ±0.87 

0.48 

1.89 ±0.67 

-0.44 ± 0.43 

030323 


RcIcJHK 

5.91 

1.43 ± 0.35 

-0.32 ±0.31 

6.35 

1.30 ±0.54 

-0.12 ±0.28 

6.41 

1.14 ±0.31 

-0.02 ±0.11 

030328 

14 

UBVRcIc 

0.23 

0.51 ±0.06 

0.00 ±0.03 

0.22 

0.49 ±0.15 

0.01 ±0.07 

0.18 

0.36 ± 0.45 

0.05 ±0.15 

030329 


UBVRcIcJH 

0.10 

0.30 ±0.22 

0.54 ±0.22 

0.09 

0.32 ±0.21 

0.50 ±0.20 

0.06 

0.41 ±0.17 

0.39 ±0.15 

030429 

15 

VRcIcJHKs 

2.60 

1.51 ±0.09 

-0.28 ±0.07 

7.72 

1.11 ±0.33 

0.05 ±0.19 

1.79 

0.22 ± 0.24 

0.40 ±0.10 

030723^ 


UBVRciJsHKs 

0.16 

0.58 ±0.26 

0.48 ± 0.29 

0.15 

0.58 ±0.25 

0.42 ±0.28 

0.19 

0.66 ±0.21 

0.32 ±0.22 

040924 


VRcIcK 

0.10 

0.59 ±0.61 

0.21 ±0.62 

0.10 

0.58 ± 0.64 

0.22 ±0.62 

0.09 

0.63 ± 0.48 

0.16 ±0.44 

041006 

16 

BVRcIcH 

0.06 

0.36 ±0.27 

0.11 ±0.23 

0.05 

0.32 ±0.33 

0.14 ±0.28 

0.05 

0.34 ±0.30 

0.12 ±0.23 


^denotes the corresponding number in the sample constructed in Paper 11. 

^Filters that are not used for the fit (e.g, due to Lyman a dampening) are not included. The degrees of freedom of the fit are always number of filters minus 
three, except for GRB 000131. 









^For GRB 980519 and XRF 030723, redshift estimates were employed (see § 2). 

‘^With only three viable colors, GRB 000131 could not be fit with three free parameters. The fit was performed without extinction. No dust preference can 
be derived. The fit has one degree of freedom. 
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Table 2. The SED Golden Sample constructed from Table 1 


CRB 

Dust^ 

X^dof 

/5 

Tlv 

970508 

M 

2.73 

0.32 ±0.15 

0.38 ±0.11 

990123 

S 

0.20 

0.46 ±0.12 

0.04 ±0.05 

991208 

M 

0.20 

0.23 ±0.37 

0.80 ±0.29 

991216 

M 

0.12 

0.38 ±0.11 

0.13 ±0.08 

000131 

— 

<0.01 

0.66 ± 0.34 

0 

000301C 

S 

1.14 

0.59 ±0.12 

0.12 ±0.06 

000911 

M 

0.26 

0.75 ±0.26 

0.20 ±0.22 

000926 

S 

0.37 

1.01 ±0.16 

0.15 ±0.07 

010222 

S 

0.58 

0.76 ±0.22 

0.14 ±0.08 

011121 

s 

0.49 

0.61 ±0.13 

0.39 ±0.14 

020405 

s 

1.75 

0.94 ±0.19 

0.15 ±0.16 

020813 

s 

0.95 

0.81 ±0.14 

0.12 ±0.07 

021004 

s 

0.69 

0.67 ±0.14 

0.14 ±0.05 

030226 

s 

0.32 

0.57 ±0.12 

0.06 ±0.06 

030328 

s 

0.18 

0.36 ± 0.45 

0.05 ±0.15 

030329 

s 

0.06 

0.41 ±0.17 

0.39 ±0.15 

030429 

s 

1.79 

0.22 ±0.24 

0.40 ±0.10 

040924 

s 

0.09 

0.63 ± 0.48 

0.16 ±0.44 

041006 

M 

0.06 

0.36 ±0.27 

0.11 ±0.23 


^This represents which dust model is used in the 
fit. Bold M and S mean a MW and SMC dust fit, 
respectively, with strong arguments favoring this fit. 
M and S mean a MW and SMC dust fit, respectively, 
with weak arguments favoring this fit. If no letter is 
given, then Ay = 0, and no conclusions can be drawn 
about the dust model. 
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Table 3: Nh along the line-of-sight to the GRBs 


GRB 

Nh (lO^^cm 

Method^ 

Reference 

970508 


X 

Galama & Wijers (2001) 

990123 

5.4ir. 

X 

Galama & Wijers (2001) 

991216 

6.8l®'3 

X 

Ballantyne et al. (2002) 

000301C 


0 

Jensen et al. (2001) 

000926 


X 

Piro et al. (2001) 

010222 

124^ 

X 

Stratta et al. (2004) 

020405 

4 7+3.7 

X 

Mirabal et al. (2003) 

021004 

< 0.1 

X 

M0ller et al. (2002) 

030226 

q 9 + 1.5 

X 

Klose et al. (2004) 

030329^ 

< 0.2 

X 

Tiengo et al. (2004) 

030429 

4.04?;® 

0 

Jakobsson et al. (2004a) 

041006 

q 9-I-0.16 
'^•^-0.16 

X 

Butler et al. (2005b) 


“Method: X - via X-ray data; O - via Lyman a (optical) 

^The value for GRB 030329 has been transformed into an upper limit. 
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Table 4: Redshift corrections of the afterglows listed in Table 2 


GRB 

dRc^ 

GRB 

dRc 

GRB 

dRc 

GRB 

dRc 

970508 

-0.16 

000301C 

-2.09 

020405 

0.77 

030329 

3.87 

990123 

-1.22 

000911 

-0.49 

020813 

-0.83 

030429 

-3.56 

991208 

-0.35 

000926 

-2.36 

021004 

-2.55 

040924 

0.12 

991216 

-0.27 

010222 

-1.32 

030226 

-1.84 

041006 

0.64 

000131 

-3.64 

011121 

2.08 

030328 

-1.09 




is the second term of Equation A5. 
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Table 5: Absolute magnitudes Mb of the afterglows 


GRB 

Group®" 

Mb (at one day'’) 

Mb (at four days° ) 

970508 

II 

-2r85lS;“ 

99 p:ro+0.14 

zZ.OO_Q 

990123 

I 

90 91 +0.03 

ZO.Zl_Q Q5 

-20.921°;°^ 

991208 



-23.251°;}® 

991216 

I 

-24.89l°;}J 

-22.721°;}} 

000131 

000301C 


—OK ni +0-09 

^O.Ul_Q Q9 

90 qpr+O.OO 

ZO.OO_Q Q9 

000911 



91 1 q+0.26 

zi.io_o 30 

000926 

I 

OK ^ 9+0.08 

Z0.1Z_Q QQ 

- 21 . 80 IH 9 

010222 

I 

90 yo+O.OS 
^ 0 . /0_Q Q7 

-21.851°;°? 

011121 

II 

-22.701°;}® 

-18.921°;}® 

020405 


-22.391°;2o 

-20.031°;}° 

020813 


90 90 +O.O 8 

^O.ZO_Q QQ 

91 1 q+0.08 

Z1.10_o 09 

021004 


-25.581°;°} 

-24.071°;°} 

030226 

II 

—22 QQ+°-°® 

zz.aa_Q QQ 

-19.541°;°? 

030328 


-21.95l°;J} 


030329 

II 

-23.931°;}® 

-22.201°;}® 

030429 


-24.921°;}^ 

-20.751°;}? 

040924 


-19-841°;}^ 


041006 


-20.761°;}^ 

-18.891°;}? 


“Denotes membership in group I or group II of Gendre & Boer (2005) 
^after the GRB, assuming at 2 :=! 
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Fig. 1.— The spectral energy distributions of all GRB afterglows presented in Table 1. 
Included are the fits to a Milky Way (MW) dust model (dotted), a Large Magellanic Cloud 
(LMC) dust model (dashed) and a Small Magellanic Cloud dust model (solid). This includes 
fits with unphysical results (e.g., negative extinction). Beneath the designation of the CRB, 
the filters are given that the SED is based on. The absolute value of the flux density is not 
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Fig. 1.— continued. 
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Fig. 1.— continued. 
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(magnitudes) 


Fig. 2.— The distribution of the derived host galaxy visual extinction Ty in the source 
frame for the bursts of the Golden Sample (Table 2). The data point at Ay = 3.8 is the 
lower limit on the visual extinction derived for GRB 970828 (Djorgovski et al. 2001). There 
are no afterglows with Ay >0.8, creating a Dark Burst Desert. The inset shows a zoom 
into the Ay < 0.8 region with smaller bins. The distribution peaks around Ay = 0.15 mag, 
showing that there is typically a small but definite amount of extinction along the line of 
sight to GRB afterglows. 
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Fig. 3.— The derived host galaxy visual extinction Ay in the source frame for the Golden 
Sample bursts (Table 2) plotted as a function of the redshift 2 ; (taken from the compilation 
of Friedman & Bloom 2005). A clear trend toward lower extinctions at higher redshifts is 
visible. The preferred dust models are indicated by symbol form (and color). Definite MW 
dust is only found for GRB 970508. For low redshifts, the distinction between the three dust 
models disappears. Extinction was fixed to Ay = 0 for GRB 000131. 
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Fig. 4.— The star formation rate (SFR) derived from the submm flux density at 850 /rm for 
selected GRB host galaxies as a function of the visual extinction Ty along the line of sight. 
There is a correlation, although the very large error bars prevent a deeper analysis. GRB 
010222 lies beyond the Icr error region of the lit and was therefore not included in the lit. 



ooom' 


SFR = 220 ± 130 + (480 ± 160) A^' 
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Fig. 5.— The distribution of the derived intrinsic spectral slopes P for the Golden Sample 
bnrsts (Table 2). A broad peak aronnd P = 0.7 is visible (implying p = 2.4 for a wind 
environment with the cooling break lying blueward of the optical), but several GRBs (GRB 
970508, GRB 991208, GRB 991216, GRB 030328, GRB 030429, GRB 041006) have p < 0.4 
(implying p < 1.8 for the case stated above). 
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(magnitudes) 

Fig. 6.— Dust-to-gas ratios in the host galaxies along the line of sight to the GRB of the 
Golden Sample (Table 2), with Nh values obtained from the literature (Table 3) and the 
Nh-j4v relations taken from Bohlin et al. (1978), Predehl &: Schmitt (1995) and Hopkins 
et al. (2004). Most lines of sight to the GRBs have an even lower dust-to-gas ratio than 
the SMG. This is especially evident for GRB 990123 for which a powerful reverse shock was 
observed (Akerlof et al. 1999). The only two lines-of-sight to the GRBs with high dust-to-gas 
ratios are those of GRB 021004 and GRB 030329 - in both cases, the Nh values are only 
upper limits. Yet both SEDs are best fit by SMG dust. Incidentally, these are also the two 
afterglows with the most pronounced fine structure (Paper II). 
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Fig. 7.— The observed i^c-band light curves of all 30 afterglows in the sample of Table 1. The 
data have been corrected for Galactic extinction and host galaxy contribution. In afterglows 
where supernova light dominates at late times, these data have not been plotted, except for 
GRB 030329, where the light curve has been corrected for the supernova contribution. At 
one day after the burst, there is a spread of 7.5 magnitudes, with the afterglow of GRB 
030329 being the brightest and the afterglow of GRB 021211 the faintest. The two inclined 
lines are meant to guide the eye, corresponding to a decay slope of a = 1.2. 
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Fig. 8.— The distribution of the apparent magnitudes of the afterglows of the Golden 
Sample (Table 2) after shifting them to a common redshift of ^ = 1 (§ 3.6). This allows a 
direct comparison of the light curve evolution and luminosity. Compared with Figure 7, the 
magnitude spread has been reduced. It is now 5.7 magnitudes at one day after the burst. 
The brightest afterglow is GRB 021004, but it is possible that the afterglow of GRB 991208 
was even brighter at earlier times when it was not yet discovered (Castro-Tirado et al. 2001). 
A large number of GRBs have afterglows of similar brightness, these have Rc ~ 19 at one 
day after the burst. Nine of the 16 afterglows that have data at one day he in a region only 
two magnitudes wide. The two inclined lines are meant to guide the eye. The vertical lines 
lie at one and four days after the burst, the times when we derive the absolute magnitude 
Mb (Figures 9, 10). For comparison, the inset shows 748 quasars taken from the third SDSS 
quasar catalog (Schneider et al. 2005) at .2 ~ 1, having r ~ 19. If GRB 990123 had been at 
z=l its optical flash would have peaked at R = 7.6, which is ~ 35, 000 times as luminous. 
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Mg (one day after the GRB at z=l) 

Fig. 9.— The distribution of the absolute magnitudes Mb, as it follows from Fig. 8 at one 
day after the GRB. The mean value is Mb = —23.2 ± 0.4. Note that only 16 afterglows 
have observational data at that time, GRB 991208, GRB 000131 and GRB 000911 are not 
included. 
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Fig. 10.— As Figure 9, but at four days after the GRB at ^ = 1 (Table 5). The mean value 
is Mb = —21.4 ± 0.4. Note that only 16 afterglows have observational data at that time, 
GRB 000131, GRB 030328 and GRB 040924 are not included. 
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Fig. 11.— The distribution of absolute magnitudes Mb at one day after the GRB at z = 1 
(Table 5), divided in GRBs with .2 < 1.4 and .2 > 1.4. The means are Mb = —22.4±0.6 and 
Mb = —24.1 ± 0.5, respectively. 
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Fig. 12.— As Figure 11, but at four days after the GRB (Table 5). The means are Mb = 
—21.2 ± 0.5 and Mb = —21.8 ± 0.6, respectively. 



























